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ABSTRACT 


This thesis explores the use of a pointed seeker lens 
Gesigned using a spherically symmetric gradient refractive 
index (GRIN). The design helps to solve the current design 
conflict between optical quality and aerodynamic drag in- 
herent in hemispherical seeker lenses. Equations for lens 
design and the evaluation of off-axis lens performance have 
peen develoved for both a homogeneous version and a GRIN 
version of the pointed seeker lens. The homogeneous lens 
1s used aS a comparison and a check for the GRIN lens. A 
FORTRAN program (GISL) has been written and employed to 
evaluate and compare both the homogeneous lens and many dif- 
ferent configurations of possible GRIN lens designs. Results 
indicate that the GRIN lens has highly superior off-axis imag- 
ing performance as compared to the homogeneous lens. The best 
results were obtained for the GRIN lens with a fifty percent, 
oositive, spherically symmetric gradient index with center of 
symmetry interior to the lens. Only very slightly inferior 
performance was observed with a five percent version of the 
same lens; such a lens possibly can be manufactured today. 
GRIN lens performance also indicates that for objects off- 
axis by more than 17.2 degrees a large scale, multiple ele- 
ment sensor array may be required; with such a sensor array, 
objects off-axis by more than 37.2 degrees may require mirror 


elements to compensate for image movement. 
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I. INTRODUCTION 


Historically, the design of tactical missiles employing 
passive or semiactive infrared (IR) seekers has involved a 
difficult compromise between aerodynamic requirements and ©) 5) 2a 
cal or seeker requirements. Whereas aerodynamically the mis- 
Sile nose region should be sharp in order to reduce drag, 
optically it should be hemispherically shaped for image 
quality and as large as possible to increase aperture and 
therefore acquisition or tracking range. Some IR homing mis- 
Siles designed for very short range anti-armor missions have 
totally ignored nose drag in order to optimize seeker per- 
formance while other designs for longer range missiles re- 
quiring high cruise velocities and greater aerodynamic 
efficiency have used the blunted ogive as a compromise. There 
have not been any IR designs which have ignored optics in favor 
of aerodynamics; nor has there been employed a pointed seeker 
lens with the desired optical qualities. 

In order to increase the performance of optically guided 
missiles beyond the current state of development, the conflict 
between aerodynamic requirements and optical restrictions must 
be resolved. Significant improvement in missile performance 
by increasing thrust is not likely due to the highly advanced 
State of propulsion today. One way to resolve the aerodynamics- 


Optics vroblem is to design a pointed lens which has, if not 
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imaging quality, enough Optical performance to allow the 
BeGuceion of tracking data. 

The lens must be in the general shape of a cone or Ogive 
with half angle compatible with the design free stream Mach 
number in order to maintain an attached shock wave. Such a 
spike-shaped lens might also be used in the diffuser pertion 
of an optically guided ramjet with nose inlet to conserve 
stagnation pressure. A diffuser-lens may benefit by a semi- 
isentropic spike shaped lens. 

Poor optical performance due to the pointed shape of 
sharp lenses has precluded their use. Gradient Refraceave 
Index (GRIN) materials, however, permit the lens designer 
the freedom to spatially vary the lens index of refraction 
to compensate for a traditionally poor optical shape. Al- 
though the use of GRIN has seen widespread use in fiber optics 
technology, it has been used infrequently in lens applications 
until recently. Lens designers are discovering that multiple- 
element photographic objectives may be redesigned using a 
two-element gradient lens [{1l]. It should be noted, however, 
that such GRIN lens have not yet been successfully fabricated 
even though large index changes in glass have been accomplished 
by the diffusion of doped electropolarizable ions [2]. Ex- 
tensive research is being conducted in the creation of ever 
larger and more precise gradients. A total change in refrac- 
tive index of approximately five percent is currently 


attainable. 


a 


At the Naval Postgraduate School, Frazier [3], Terrell 
[4], and Amichai [5] have studied GRIN as applied to the 
sharp lens problem. The brief introduction by Frazier was 
followed by Terrell who designed a sharp lens having a coni- 
cal outside surface and a variable inner surface using a 
homogeneous index (HIN) and then also briefly touched on the 
GRIN application. Amichai extended Terrell's lens to the 
GRIN case and included preliminary results from this thesis 
in a computer routine intended as a framework for lens 
optimization by following researchers. 

This thesis investigates the design and off-axis per- 
formance (skew rays) of a GRIN seeker lens having a variable 
outer surface and a fixed right circular cone as an inside 
surface. The variable outside surface is determined by the 
character of the spherical gradient employed and varies from 
a pseudo ogive to an approximate isentropic spike as shown in 
Figure l. First, the theory of the HIN lens is developed, 
followed by the design of a homogeneous lens intended both 
as a comparison and a check for the GRIN lens theory and the 
design which follows. Lens performance parameters are dis- 


cussed and results presented for both the HIN and GRIN lenses. 
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OUTSIDE CONE 


faye Geeretl’s Lens With Conical Outer Surface 


mio DE “CONE 


(b) Lens Shapes In This Thesis With Conical 
Inner Surface. 


Figure 1. Lens Shapes 


ae, 





Il. THE HOMOGENEOUS LENS 


=  DHEORY 
Snell's law is the cornerstone of contemporary lens de- 
Sign in homogeneous optical materials. In the HIN case 


Snell's law is used in the familiar form 


N. sin @ = alg 8 Cis) 


where N. Homete index Of refraction of the material corres- 
ponding to surface intercept angle Os with respect to the 
Surface normal, at the interface between surfaces. Lens 
geometry and the relative values of N; determine the resultant 


optical behavior. 


B. ASSUMPTIONS AND SIGN CONVENTION 

In order to simplify the design and analysis of the 
seeker lens problem, certain assumptions have been intro- 
Guced. Although energy loss upon transmission through the 
lens at each surface is calculated, it 1s assumed that the 
light is monochromatic radiation, time dependent 
electric and magnetic fields. Light impinging upon the lens 
is assumed to have a planar wave front as if propagating from 
an object at infinity. The presence of a shock wave attached 
to the lens is ignored as are any other regions of expansion 
Ome Ceampression in the flow field about the lens [6]. Further- 


more, the index of refraction of the free stream 1s assumed 
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to be equal to that of the interior space behind the lens 
even though computer routines were written with the flexi- 
bility to process unequal values. 

The sign convention used is a right handed system with 
Spatial coordinates positive to the right, up, and out of 
the page as seen by the reader. All angles are assumed 
positive counterclockwise from point of reference unless 
otherwise noted. All linear dimensions are implicitly non- 
dimensionalized with respect to the maximum radius of the 


tagmt Ciucular cone forming the inside surface of the lens. 


= wet LENS DESIGN 

i= lens Gesign procedure consists of calculating both 
the loci of points forming the outside surface and the slope 
©f the surface at each point in the meridian plane. It is 
convenient to approach the problem by placing a point source 
of light at the design focal point F and calculating succes- 
Sive refracted ray paths U' using Snell's law at points Ds 
on the inside surface in Figure 2. Points Ee are formed by 
tne Intercept of the refracted ray and the slope of the out- 
side surface as extended from the previous point =,_,. The 
tangent to the surface which will refract the ray in the de- 
simed direction, U, may again be found from Snell's law. 
Gee weecuracy Of the calculated surface increases as the 
spacing between points decreases or as angle AU" becomes very 
small. The final ray is parallel to the lens axis and point 


a Gemnecaponds to point B, which is the origin. The area 
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Gommed by the triangle ABE | is Opaque to prevent rays from 
Beaetering through the opposite side of the lens. 

® Parercular lens is begun by specifying the focal 
length BF, the cone half-angle a, the indices of refraction 
Ny: No: and N3, the cone radius at the edge R, the ray 
direction angle U, and the total number of rays to be traced. 
The thickness, T, of the lens at the edge (E,D,) must also 
be specified in order to define which of the family of possi- 
ble outer surfaces will be calculated. In general, angle U 
will be taken to be zero in order to investigate objects on- 
axis at infinity. R will always be set equal to one and to 
nondimensionalize all linear dimensions every length is im- 
Pbueicdely expressed as a ratio with respect to R. The angle 
8 at the focal point is measured from the lens axis to Dy and 


ls expressed as 
tan |8| = qgpoacotan ay 
Gimges if D+) is the total number of rays to be traced, then 
SS = Get Sy, 
and 
Uy = 6 = J-Av" (4) 


where J is the ray number. To begin the lens design the 
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coordiantes X(T) vy, (J) eng ele abele Ds must be found. Here 
the subscript 2 refers to the inside surface where x, (9), 
y;{9) ecm COordinates of E; On the outside surface. 


Following Kingslake [7], Q in Figure 3 may be expressed 


as 
& = BF sin U"“ = x5 (J) sin U" + Y5(J) cos U" (Ss) 
Buc, 
es) = = x,(J) tan a (6) 
so that 
BF sin U" = x5 (J) [sin Oe — tan a "cos Ur] (7) 


Now, x5 (J) and Y5(d) are 


2 BF sin U" 

Xo (5) Sets tama cos UU") (8) 
= Been U" tan a 

YJ) lo (Sener ean a cos u | (9) 


EeokKing now to find the incident angle I, which the ray forms 


with the inside surface at Ds Lt can be seen that 


BF - x5 (J) 


tan (a + I5) = years o (10) 
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Substituting into Equation (1), Snell's law, the refraction 
angle I. can now be written as 


Sie 
ya Uae) sin 15] (12) 


The path of the ray inside the lens may now be expressed 


in the form of the angle U' where 


- 


5 (13) 


To this stage all rays are treated the same. Now, however, 
a differentiation must be made between the outermost ray 
which defines the edge of the lens and all other subsequent 
rays. 

For the outermost ray the thickness, T, of the lens at 
the edge can be specified which yields x, (1) and y, (2) 


immediately: 


ota x, (1) = CoS. U' (ra 


y, (1) y5 (1) + T cos U' (15) 


At point E, we may write Equation (1) again to obtain 


il 


fae ance 1a iors (16) 
IT, = sin No Sa. il 
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where angles Ty and Ty are the incident and refracted angles 
at the outside surface and are not yet known. Angles I 
and Ti may be found by first noting that 


U+ tI = U' + Ty (ale) 


then substituting Equation (16) into (17) and rearranging 


temms to obtain 


N 


1 jee 
-_" I,J (18) 


= ' - see 
I, U Ue+ Sin 


It can be shown after some algebraic manipulation that 


Equation (18) may be solved for the angle I, inehes forum 


= See cue ein) a 


[cos (U'-U) ~ N)/N,] + Sin (U' — U) 


eegre i- may be found, if desired, by substitution of Equa- 
tion (19) into Equation (16). Of more importance, however, 
is the determination of the slope of the outer surface and 
the surface normal. Now that angle I is known these slopes 
may be written as 


dy 
dx!tnormal 


DYDXN (J) - tan(U + I,) (20) 
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: - dy _ 
bY PRT (J) aolccntoce.. = cotan(U+L 


wange»ne 


1? C2) 
Equations (20) and (21) may be evaluated by using Equation 
£ 
Co) Bor I,- 
Now that all the parameters are known for the first ray, 


the remainder of the points E_ and the respective slopes 


J 
may be found. Each successive ray is traced as before by 
Equations (4) through (13). Equations (14) and (15), how-~ 
ever, may not now be used since the lens thickness along the 
ray can not be specified. Instead, the intersection of the 
ray and the surface slope from the previous ray 1s used to 


define the new point E The intersection is found by first 


“Je a 


writing the equations of lines representing the ray and the 


surface tangent. For the ray: 


w= -x% tan U' + Y5 (J) + x5 (J) = ag) 2) (22) 


mor the surtace tangent: 


y = xX cot (I, + U) + y, (J-1) - x, (J-1) cot(I, +) 
Cs 


Equations (22) and (23) are solved simultaneously to yield 


the coordinates of E. which are 
W 





y,(9) = { } (24) 
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where: 


A = cot(I,+U) Seen) Ly) eo) amcor) | (25) 

B = ya) = x, (J-1) cot (I, + U) (26) 

f= 1] + cot(I, + U) COmCu™) (277) 
and 

x, (J) = cot(U") [-y, (J) ty, (J) +x, (J) tan(U') | (28) 
ow Equations (19), (20), and (21) are used to calculate I 


1 


and the slopes at E Therefore, the remainder of the lens 


3° 
surface may be generated. The opaque region at the surface 
1s formed by extending the slope of the surface at Ex to 
intersect the lens axis. Here, K refers to the last ray. 
The nose half-angle, y, thus formed is given by 


tany = & (29) 


surface, Ey 


D. SKEW RAYS 

A skew ray is one that begins from an off-axis object 
point and enters the lens either in front of or behind the 
meridian plane (z = 0). For every skew ray, there is a 
corresponding mirror image skew ray on the opposite side of 


the meridian plane so that two skew rays are traced at the 
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expense of only one calculation. These two skew rays inter- 
Mece ate the same diapoint. 

Large numbers of skew rays are traced through a lens 
in order to study lens performance at different obliquities. 
The procedure is to superimpose a grid over the lens aperture 
and to trace rays through the intersections of the grid, 
through the lens and onto the image plane. An image plane 
spot diagram and an energy density plot are then constructed 
fer study. 

To accomplish this, the aperture grid has been attached 
to the nose (opaque region) of the seeker lens at station A; 
see Figure 3. The plane of the aperture grid is tilted 
relative to the lens by a variable angle. A transformation 
between the grid coordinates and the lens coordinates has 
been derived to connect skew rays from grid to intercept with 
the outside surface of the lens. Referring to Figure 4, it 


can easily be seen that 


x" = x + ABIGoS a_- sin a (30) 
; oe p 

= + AB|]sin toy COS i (32>) 
V [x Jsi on \ . 
ie =) Z (32) 


where re ie the tilt of the grid plane, AB is the length of 
the opaque nose portion on the x-axis. Both z and 2’ are 


Pes tive out of the page. 


20) 
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The three dimensional lens outer surface is generated by 
rotating the array of outer surface coordinates in the merid- 
maeplame through 2n about the lens axis. Each pair of 
coordinates x (J) ,y, (9) thus describe a circle in the y~z 
plane of the lens. This circle transforms, however, to an 


ellipse in the grid plane given by the equation 
are) x + oeB)}) sin O65 z VR? - z'* cos Bo (3) 


where R is the radius of the circle in the x-y plane. R 


may be expressed by the familiar equation of a circle 
R == ge eZ (34) 


in the y-z plane of the lens, or by solving Equation 3.3) 
mere RR in the grid plane 


2 1G) 


R = z'* + (eee ~ (x + AB) tan oe: 
With the aid of the foregoing groundwork, the x,y,2 coor- 
dinates of the ray intercept with the outer surface may be 
found. To see how this is accomplished, first refer to Figure 
5. The skew ray will pass outside the circle formed by rotat- 
ing some point x, (J) +y, (9) and inside the next circle formed 


by petating x, (J-1) ,y, (3-1). In so doing the skew ray will 
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intercept an imaginary cone formed by the two emolecn. 


Now, the y',z' of the skew ray in the grid plane and a 

trial value x, (J) from the array of surface points are sub- 
stituted into Equation (35) and a value for R is calculated. 
If Ris greater than Y; (9) but less than y, (3-1) than the 
Sereeputace Circles have been found. If not, another trial 
value of x, (J) 1s picked. Inpractice the x4 (J) corresponding 


to 
y, (9) a (36) 


is chosen as the first trial value and subsequent trial 
values are x, (J-1), x, (J-2) anGesoeroren.. 

The equation of the cone passing through the two circles 
may be derived by again referring to Figure 5 and noting 


that 
ep) Sop!’ = Vy 4 2° = R (S70 
OP = y = m(x - x, (J) J 1p yd) 37) 


where the slope m is given by 


trecall that a straight line between points E, was used 
to form the outer surface in the meridian plane. These 
approximations are valid only if the spacing between points 
or circles is very small. 





y, (J-1) ~ Y; (J) 


ee x, (J=1) ~ x, (J) a 


and therefore 


RU = UU (aaty =, cay! EX XG) iF y,(g)} (2s), 


by C@uatang Equation (39) to Equation (35) and solving for 
x, the expression for the x-coordinate, Xo of the outside 
surface ray intercept may be found. After some algebra, this 


reduces to the complicated relationship 


PARL , PARI _ PAR3 





") ~~ BARD ~  $$PAR2 £PAR2 (40) 

where: 
2 Z 

PARL = ab+ce-b AB-cd (4) 

Dame = BS = (42) 

pe? = 62 b> (AB) - = 2abaAB + z' - (e=cd) - (43) 
and where 

a = y'/¢os Bs (44) 

> = tan a (45) 





Y, (5-1) -y, (J) 


a ig WHT) =x, (46) 


Qu 
NH 


x, (d) a7) 


m 
I! 


Y, (9) (48) 


It has been found that the plus sign in Equation (40) gives 
the correct values until PARL/PAR2 becomes greater than 

x, (1) at which time the negative sign must be used. It now 
follows from Equations (31), (32) and (40) that the Yo and 


Zo coordinates of the intercept are 


Sa 
Yo Boe a (X) + AB) tan os (49) 


and 


Z = Zz' (50) 


The next step in tracing the skew ray is to ascertain the 
direction cosines of the ray inside the lens after refrac- 
tion at the outside surface. This is accomplished in a 
series of steps beginning with the determination of the 
angle ¢ which the skew ray makes with the normal to the 
Sltemece at the point of intercept. Angle $¢ may be found by 


meemnGg the scalar product of the direction cosines of the 
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skew ray and the surface normal. Thus to find od, the surface 
normal must first be obtained. 
As might be surmised by the reader, the direction cosines 

of the surface normal can be found by taking the gradient 
of Equation (39) with Equation (37) substituted for R. 
several orders of magnitude may be gained in accuracy, how- 
ever, if the normal Na is found by interpolation in the 
meridian plane and then rotated to the point of intercept by 
angle w as in Figure 5. This interpolation procedure involves 
multiplying the ratio S/ST with the difference between 
DYDXN(J) and DYDXN(J-1) and adding the product to DYDXN(J). 

DYDXN, 4 = S(DYDXN(J-1) - DYDXN(J)) + D¥DXN(J) (51) 
were subscript PIM refers to point of intersection, meridian. 
Values of the slopes of the normals in the meridian plane, 
BYDXN(J), are given by Equation (20) and § and ST are the 
linear separation of points P", PIM and P", P respectively. 
myo. 1s given by 


5 ) eZ 
— I(R- y, (J)) - One a ] ae 
ST (anse . tech? c 


Mherefore, the normal vector, Nee, in the meridian plane is 


a 


= aoe) 1 t+ C”7F (555 
m DYDAN Dow 





and since 


tan yp = =— = — (54) 


the normal vector at the point of intersection on the 


outside surface becomes 


N = (Sypma + +eGos( tan i le + sin(tan Bie (S53) 
PIM YQ YQ 
and 
_ 25 WD 
iIN'| = (VAP ee (36) 


as 


Now the scalar product of the skew ray vector, R, and surface 
normal, N', may be performed to find ¢. Since the acute 
angle between these vectors is required, the dot product 


must be written 


> 


R°N' = |N'! cos(1-9) (57) 


| 


Substituting Equation (55) for N' and noting that 


“~ ~ 


R = Ki+Lj)+Mk = cos ani - sin a, (58) 


mae beet side of Equation (56) may be expanded to 


a ces a lh a (59) 
R°N' fee Sin cg Ccos(tan =) 
DYDAN Dy p Yo 
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Solving Equation (56) for the angle ¢ land sineroducs ng 


Equation (58), the expression for angle 6 becomes 


COS @ 7 
DYDXN - Sin a_ cos(tan ——) 
“ roeog"? | ——_"pr_? iy _| 
3 Z _ ee. | 
[DYDXN,,_., 10) 4 sing? (tan ) —9) ) 1/2 
Yo Yo 


Z 


fee (CNA 


(60) 


Now that the acute angle, $, between the ray and the outside 
surface normal is known, 9$', the angle between the refracted 
ray inside the lens and the surface normal may be found using 
Snell's law. 


N 


* (= siya oh (61) 
2 


\' = “sane 
The next step is to ascertain the direction cosines K', L', 
M" of the refracted ray inside the lens. Following Kingslake 
[7], Figure 6 shows the optical vector relationship between 
che skew ray, R', the refracted skew ray, N', the surface 
normal, and the indices of refraction of the two media. 
Algebraically this relationship is written as 
cos )N' (62) 


N.R' = N 


5 12 + (N, cos > -N 


i 


By sesolving Equation (61) into component form, the direction 
cosines of R' may be found. After rearranging, K', L', and 


M' become 
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No 


N, 





Flqure 64 Kingslake’s Skew Ray Diffraction 
Diagram 
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N N 


i 1) 
K - n, X + (cos 9' - 7 cos $)k (6.3) 
Ny Ny 
Ii = a elealices ¢'= y COS o)2% (64) 
2 2 
Ny Ny 
mM” = ~=-M+ (cos 9' - = cos $)m (65) 
a, ay. 


where K, L, and M are given by Equation (57) and k, &%, m are 


-1 
_ _ BYDENery 
ae (66) 
[Ne | 


costeaan — —-) 
y 
L = (67) 


(68) 


Again, |N in Equations (66) through (68) is expressed by 
Pe@ation (56). 

At this stage in the process of tracing the skew ray, 
the coordinates of the external ray intercept Xgr Yor 2 
are known as are the direction cosines K', L', M' of the 
skew ray inside the lens. The next step is, of course, to 
find the ray intercept with the inside surface (cone) and 
Bmew dissection of the ray subsequent to refraction. The 


intermediate steps are similar to, if not identical with, 
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the foregoing. One major difference, however, is that the 
inside conical surface may be expressed analytically and the 


-, 2, May be found exactly without 


Meeercept coordinates x., y 
1 1 a 


approximation. Proceeding, it can be seen that if > re es ae 
1 a 
z, were already known, Equations (63), (64), and (65) could 


be rewritten as 


a a 0 
K as cee (69) 
age! 
eae. ae 0 
L + Sane (70) 
Z.- 2 
vi +H (aie) 


where D is the distance between surface intercept points. 


Upen rearranging: 


ee OSE Fae I, 
y; = DL' + yo (73) 
zZ, = DM' + Z (74) 


Furthermore, there exists a relationship between Xin Yas and 
Zs which is given by the expression for the inside conical 


surface: 


ye + 27 - x? tan*a = 0 (75) 
a8 ae 





where a is the half-angle of the cone. Substituting Equa- 
cies (72), (73), and (74) into Equation (75) and solving 


mem D, 1t is seen that 





ye 
=-p P, = P.p 
Mh a 
Clg aes a : (76) 
2 Z 
where 
P = Vay SMe. =  K'y Mya (72) 
ate a 8 at) 0 
Ps = g'2 4 m2 _ RO Bon g (78) 
and 
_ 2 2 2 2 
Pz = Yo + 29 - Xj tana (CaS) 


Correct values for D are obtained by using the minus sign in 
Pemaecion (76). New that D is known, values for K', L', M', 
Kye Yor and Zz are SUMSiGLCUGEaA into Equations (72), (73), and 
(74) to yield the coordinates of the inside surface intercept 
peant. 

In the case of the inside surface, the gradient may be 
used to obtain the surface normal at the point of intercept 


Since this surface has been expressed analytically. If 


Mem@ation (75) is denoted by f(x,y,2) then 


= —— (80) 
Wee 





where N; is the unit vector in the direction of the surface 


Hommal at the point of intercept. 


ees ee, 8 
VF = yx i ts ay 7 ¢ 5 
Smeal cernately, 
Ties i tof 5 + £_k 
Xx y Z 


Thus, the direction cosines of 


be expressed as 


gio= 


where: 








Day Ae 
be 
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Sen 
”~ > 


tne Normale’ | 


Here, 


(81) 


(82) 


Lam May 


(83) 


(84) 


(85) 


(86) 


Gam) 


(88) 


and 


ths Dlg s 72 
2 
ve 2 (xf tan’ a + yy - za (89) 


Now that the two unit vectors R' and N' are known the inci- 
dent acute angle of the skew ray with the inside surface, 


¢., May be ascertained by again taking the scalar product as 


ae 


Ma Bauation (57). Thus, 


o, = m=cos + [K'k' en ee ee (90) 


and using Equation (1) again 


-l No 
t —_ Q _f& ° 
?, = sin seaman COI) 
where o 1s the acute angle between the skew ray and the 
surface normal after refraction at the inside surface. 
Direction cosines K", L", M" of the skew ray after refraction 


are feund analogously to K', L‘, and M' in Equations (63), 


ha, > Mm; ois oi, 





(64) and (65). Here, however, Nos N3, 
ema k', 2', m’ are substituted for Nyy No. Ky bo i oo 
ame x, £, m respectively. Thus 
N. N. 
K" = — K' + (cos o! - — cos 9.)k' (92) 
N 1 N aE 
3 3 
N. N. 
iG = N, Pies cos do, - Ny cos >,)2 (93) 
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and 


N N 


A = t Zz t 
M = N, M' + (cos oH ~ N, cos d;)m (24) 





Finally, the skew ray intercept with the image plane may 
be found. The image plane is treated as another surface 
along the path of the skew ray, and the intercept coordinates 
ne 


im’ Yim’ Lie age Casy tow rma. Equations (72), (7/73) and 


(74) may be used again in the form 





= L vf 
Xen D'K" + x. (35) 
Yim — OF) 4 Va (96) 
Ae = D'M" + Zz. (917 ) 
im al 
where D', in this instance, is the linear separation between 
coordinates 5 vey a at the inside surface and Apa 0 ene 
a of the image vlane. Furthermore, 
x. = BF (98) 
slay 
is the equation of the image plane. It follows that 
(BF - x.) 
—— Se 
D a (22%) 
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after substitution of Equation (98) into EBQUaE1On (9 o)) % 


Therefore, 
BF = Xs 
= = jae 1 ye) (100) 
and 
BF - Xs 
os = [eer aes | a at (101) 


Hence, the skew ray has been traced onto the image plane 
meptme. focal point of the lens. The ye coordinate of the 
corresponding mirror image skew ray is the same as Equation 
(100). The Biter coordinate, however, is the negative of 
Equation (101) since the mirror image ray is behind (when 
viewed along the z-axis) the meridian plane. After tracing 
a complete set of skew rays through the lens a spot diagram 
may be plotted. Clearly, the number of rays to be traced 
depends entirely upon the incremental size of the aperture 
feria enosen. According to Kingslake [7], at least 100 rays 
must be traced to give a fair approximation of the actual 
image. In addition to the spot diagram, an energy density 
plot may now be constructed by counting the number of rays 
within progressively larger radii from the image centroid and 
Geemeplotting the mumber of rays as a function of radius. 
Here, each ray is assumed to contain a unit, nondimensional, 


an amount of radiant energy for convenience. 
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E. RADIANT ENERGY LOSS 


In reality, each ray loses intensity upon transmission 
at each interface. Of the total amount of energy contained 
in each ray,a fraction =» Will be transmitted, a fraction I, 
will be reflected, and a fraction Tn will be absorbed by the 


medium into which the ray is propagating. Since it has been 


assumed that absorption is negligible, it must be true that 
ie tle t=. wh (102) 


Furthermore, the relative amounts of transmitted and reflected 
electromagnetic energy may be calculated by the well-known 


Fresnel Equations which state the dependency of In and In 


upon the angle of incidence and the indices of refraction 


at the interface. That In and Ip 


the orientation of the electric vector with respect to the 


are further dependent upon 


geometry of ray incidence is fundamental to the boundary 
conditions which govern the form of the Fresnel relations as 
Gerived in Hecht-Zajac [8]. Since this thesis examines lens 
response to monochromatic radiation, the indices of refraction 
are not considered as a function of wavelength; further, the 
electric vector orientation is assumed to be rapidly and 
randomly changing with time. By time averaging field com- 


ponents,it may be seen that the reflectance is 


ee 2 
eT se (103) 
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where 


ane Reem meee ia (104) 


2 OO IRE 
cos eee Nas Sin’ 9) 
and 
2 7 2 ee ive 
Nei cos ¢o (Nes - sin 95) 





2 2 weee (105) 
Nes cos 9 + (Ny eat) ye 


Here, NL. 1S the ratio of the index of refraction of the 
transmission side of the interface to the index of refrac- 
tion of the incident side. From Equation (102) it now follows 


that the transmittance through the interface is 
im, = Lo= 1 (106) 


The total transmittance through the lens is simply the product 
Sr I, at the outside surface with that of the inside surface 
where 25 VS yembstacuted for 9 in Equations (104) and (105). 


Total internal reflection of the ray may occur at the 


inside surface if the incident angle becomes too large. 


FOllowing Reference (7), this occurs when angle 5 is equal 
Me er create; than 7r/2. Thus, Snell's law becomes 
ae 
sino. = = (b07) 
1 No 


ana any ray with Oo; equal to or greater than this will be 
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Peery internally reflected. In this thesis such rays are 
Pabeled “failed rays" since they fail to intersect the image 
rlane. If N, us) 8.25 and N, is 1.0, the incident angle for 
Seealeinternal reflection is 41.81°, or greater. 


=. OPTECAL PATH LENGTH (OPL) 

The optical path length of a skew ray is an analytical 
tool with which the researcher may ascertain the phase of a 
ray at the end of the path. By so doing, the image diffrac- 
tion pattern may be constructed which shows the addition or 
subtraction of amplitude depending upon relative phase. Since 
each ray must begin with the same phase, monochromatic radia- 
Hi0n is used for diffraction experiments. OPL is included 
here only as a matter of interest. The calculation of optical 
path length is simply the sum of the geometrical path segments 
@e a rey multiplied by the corresponding index of refraction 


Seeene mMearum for that segment. 
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III. THE GRIN LENS 


me  PBORY 


This thesis assumes a spherically symmetric, inhomogene- 
@us, 1Sentropic medium in which the refractive index varies 
from point to point but is independent of direction at each 
point. The refractive index is a function of the coordinates 
of the points of the region being considered. The problem 
of describing the resulting curved paths of rays in such 
GRIN materials has been solved long ago in the form of a 
Single second order vectorial differential equation. Marchand 
[9] has shown that the solution to the differential equation 
in the case of spherical gradients can be written in polar 


coordinates in the plane of the ray as 


dr 


viner’ -e2]l/2 


(108) 
T0 


Here Xo and 369 are values of r and 9 at a convenient reference 


Peasant On the ray and 'e| is a scalar constant along the ray 


given by 


Ce = ees sin Vp luoye)) 


Referring to Figure 7, r is measured from the center on 


symmetry of the index function; angles 9 and 9) are measured 


ge 









* A 
Pare 


le 





Figure 7. Grin Lens Geometry 





counterclockwise from the lens axis; angles wp and Wy are 
measured counterclockwise from the radius vector to the 


instantaneous ray direction vector. ¢ is the sign function 


expressed as 


€ = sgn Vo =e) GeO ) 


where 


l¥gl 2 1/2 (111) 


tmpebawewation (109), Ny refers to the index function n evaluated 
at Lo: The form of n has been chosen as 


ROSS n(O_/r) (lee) 


where e. Genotes the position of the center of symmetry on 
the lens axis (shown in Figure 7) and r is the radial coor- 
dinate from om to the point in question. This form allows 
the study of the effects on lens performance as 0. is) changed. 
More specifically, in order to allow the analytical 
integration of Equation (108), the expression used for 


n(O_,r) is a generalized version of that employed by Luneburg 


as described by Marchand [9]. Here 
i - Ue 2 
n = [a+ b(=-) J (iS) 
Z 


Fe 


Wwmereas Luneburg used specific values for a and b. Generaliz- 
ing the index function enables the strength, or percent change 
of the gradient to be varied as well as the algebraic sign 
of the gradient. If the parameter b is negative, a decreasing 
parabolic gradient results; conversely a positive b yields 
an increasing parabolic gradient. By visualizing a plane 
wave front passing through a spherical gradient, it may be 
seen that the negative gradient results in light rays bending 
forward at the center, as shown in Figure 7, whereas a posi- 
tive gradient has the opposite effect. This principle was 
used by Wood, as related by Marchand [9], in constructing 
Simple lenses having plane faces and a radial index. The 
Wood lens acted as a converging or diverging lens depending 
upon the sign of the gradient used. 

Pog@utaonally, it must be noted that if the parameter b 


mee equal to zero, Equation (113) reduces to 
oS Ve (114) 


which is a HIN lens having constant index of refraction. 
This fact has facilitated the correlation of GRIN and HIN 
computer trace algorithms. 

Equation (108) may be integrated by a change of variable 


using the relation 


vee (a (115) 
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which leads to the solution 


© 
do] m 


7 272 En [ne 
ys ——— - si 
a2 
a“ + 4be*/Re Va* + 4be*/Ré 


(116) 


Equation (116) gives @ as a function of r. This equation 


@em be GCasily solved for r as a function of @ in the form 


J22\e 


Zea? 
2 i 2e /r.--a 7 2 
at a* ae sin [-26 (8-89) +8in | | 


ay, Va? +4be*/R 





rels7) 


The instantaneous direction of the ray at any point r,0 
may be ascertained by using the invariance of e. Hence, 
e = enrsiny = en 


ry sin, (118) 


GeO 


and 


ne Yr. Sid Yo 


wees 0.0 
y = sin “l ae (119) 


Furthermore, the orientation of the plane of the ray may be 
easily deduced since every ray ina spherical medium is a 


plane curve lying in a plane through the center of symmetry. 


TS 


Using this fact, Marchand [9] has shown that a suitable 
conversion from coordinates r and 6 in the Dlete Onrtne nay 


to global Cartesian coordinates may be written in the form 


x = x6 4, ) 120 
=, Py (120) 

y = xr(6 “0 + 7 ) Gl 
EA do ) 

Z = x(6 a? Q ) C22) 
Xp 0 


Here Por Ug? a are the initial directs@mwW cosines of the 


ray at Lor Sgi Xo Yor and Zo ape theseameeouamCOOcamnates 


corresponding to ry 18° The parameters 6 and n are given 
by Marchand as 
n = sin §/sin Vo ‘ales. 
6 = cos6-y)n cos Vp eZ 


It should be noted that Equation (117) may become singu- 
lar for certain rays where w or Vo become very close to zero 
or t. This singularity may be more easily seen in Equation 
(123) where n becomes indeterminate as both 9 and Vo approach 
zero and/or m. In practice these conditions occur when om 
is located either far out in object space, coincident with 


B, Or on the image side of the lens; see Figure 8. Positioning 
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Vz 
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(-) OB 


+A N/N 
| Vj 
EXCLUDED “YY 
REGION DUE TO 
SINGULARITY IN 
EQUATION (108) 
O 
; Yj 
-AN/N 


Bxeclucded Recions for Center of 
Symmetry Due to Singularities 


TS 


B 
) 





the center of symmetry at these locations has therefore 


been avoided. 


B. ASSUMPTIONS 

The aforementioned assumptions for the HIN case also 
apply here. Additionally, it must be assumed, as depicted 
in Figure 7, that the GRIN lens could or soon may be fabri- 
cated from a sphere of dielectric material with the required 


spherically symmetric parabolic gradient. 


C. GRIN LENS DESIGN PARAMETERS 

In the HIN lens, the available design parameters are 
Basic. These include: F, R, T, a, U, and N.- Parameters 
available for varying the design of the GRIN lens, however, 
include those of the HIN case but expand the index of refrac- 
£10N variable No siahere Oo, a, +b, and -b. These additional 


lens design parameters greatly expand the lens designer's 


power to bend radiant energy to his will. 


De GRIN} LBNS DESIGN 

The GRIN lens design procedure, although paralleling 
that of the HIN case, is somewhat more complicated in that 
the rays are now curved and the index of refraction varies. 
Accordingly, the same design process is used but with more 
intermediate calculations required. 

The additional calculations arise since the intercept 
of the GRIN ray with the surface tangent cannot be solved 


in closed form, and an iterative solution must be used. 
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The reader will note that Equation (117) will Viera mr Le 
(9-59) is known. Hence the iterative procedure is to 
pam@ess (0-6,) based on th 

g 0 ce AtN ceorainates Xia Yip which 


are calculated as in the homogeneous lens. Thus 
ae ee 84 C25) 


in Figure 9 is used in Equation (117) to find point p. The 
tangent to the ray path at point p is now extended to inter- 
cept the surface tangent again using the homogeneous inter- 


cept relations to find x The prime superscripts 


lH’ Yin- 
indicate successive iteration values. The distance d is 
Employed as a measure of the error of point p. If dis not 
within an acceptable margin then oon is calculated based on 
Xiu/Yig and the procedure repeated to find p' in Figure 9. 
If d' is not within allowable error then the iteration con- 
Eemwes until it is acceptable. In practice, this iteration 
procedure has proved to be extremely rapid, rarely requiring 
more than three iterative steps before converging. Slight 
modifications, however, must be introduced to handle a posi- 
tive gradient. Furthermore, if 0. 1s located outside the 


lens proper, a decreasing angular increment must be subtracted 


inter- 


from oo to ensure that each radial vector Lays Lays a 


Sects the ray during iteration. These modifications to the 
intercept iteration procedure are recorded in the program 
Meeeiac £05 program GISL (for Gradient Index Seeker Lens). 


Refer to Appendices A, B, and C for a full description of GISL. 


Fes, 
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Figure 10. Expanded View, Intercept Geometry 
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With the overall iterative intercept procedure now being 
clear, the mathematical details of the GRIN lens design 
follow. Referring again to Figure 9, the GRIN ray is traced 
as in the HIN design procedure to obtain the coordinates 
x5 (FT) r¥5 (I) on the inside surface in the meridian plane. Now, 
before Snell's law can be used, the index of refraction must 


be found at these coordinates. Using Equation (113) 


Ip 2 Sy 7 
a = [a + b (==) ] cl25) 
Z 
where 
= By see (ae (126) 
io (x, (J) + OB) Y5 


and where 


R, = V(x,(1) + ie) © ee 0) 7 (127) 


+i - wil aa 1236 
Vo 1 U 36 ( ) 


ii 
U' is identical to the HIN case’ and 


luere Equation (125) must be substituted for N. in 
Mematcion (12). 
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Y>(J) 


-1 
8 = oA ae i ES 
oan LYS J) + OB 


J (Gie29) 
In Equations (126) through (129) OB is the line segment from 
0. @o 6B amd is defined as a positive quantity to the left 
of B for algebraic ease of manipulation. By substitution 
into Equation (109) the scalar constant, e, of the ray can 
now be found. Furthermore, since angle U' is now known 
coordinates XiurYiy may be calculated using Equations (14) 
and (15) or Equations (24) through (28). 

Por the first ray defining the edge of the lens, Equa- 
tions (14) and (15) are employed, and, since no iteration 
= required, angle (6,799) is used in Equation (117) to find 


r immediately. Here 


x + OB 
aoe a (130) 
A 1s 
‘ai 
and 
| eee ee) (131) 
H 1H 1H 
Thus 
= aoe ie 
x, (1) r cos 9, OB ( ) 
and 


oe) 


y, (1) = xrsin 8 mn Cakes) 


are the coordinates of the first point on the outside 
surface. Since e is known, Equation (119) is employed to 
find the angle ¥. Angle ¢ in Figure 10 is employed to 
translate y into the ray direction with respect to the lens 


axils by the relation 


a te a on) (134) 


meeretore, by substituting ¢ for U' in Equations (19), (20), 
and (21), DYDXN(1) and DYDXT(1) may be found. Here again, N, 
in Equation (19) must be replaced by n5, as given by Equation 
(113) evaluated at x, (1) ,y, (1). 

Now that the first point on the outside surface is known, 
along with the surface tangent, the remainder of the K number 
of rays may be processed to yield the balance of the outside 
Surface. Each subsequent intercept between ray and tangent 
must be iterated. Thus, unlike the first ray, once Xyu7 You 
Cue ie, ¥, amd Cc are Known, Xu and Ying. are ascertained by 
the substitution of the coordinates of point p for x5 (5) ry (5) 
Me Semee2ons (24) throush (28) with U' replaced by ¢. The 


Seercinates of point p are 


* 
II 


r cos é, ~- OB Res) 
Pp Jt 


Y = cS Sin 8 (136) 


'O 
ae 
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Therefore 





a C (case) 
where 
A = cot(T, + U) G6et £ Ss + x, tan ¢] (138 
B = y, (3-1) = x, (J-1) cot (I, + U] (26) 
Cc = 1+ cot(I, + U) eeites (139) 
and 
2h ees) 
i = Sin SEEN nEESSEREDEEESEREEY SAID SPE (140) 
Lcos(¢-U) - N,/n,] + sin’ (T-U) 


Note tnat n. in Equation (140) is found from Equation (113) 
evaluated at ro: 
Additionally, 


Xu = cot 6 E-Yin + his + aS tan 7] (141) 


Now the error, d, may be evaluated as 


ee 1 ce 2 
d= /(x) > yy + yy! 7 Yo) (142) 
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Here d 1s compared to ieeeto >: If d is larger than this 
value then the entire procedure is repeated by substituting 


Xe Yin ' for Xu Yig and so forth. Once the error criteria 


are satisfied 


: 
va 


x, (J) Oe) 


y, ‘9) (144) 


iI 
he 


and the next ray is processed. Correlation between GRIN 
(with b set equal to zero) and HIN design procedures run 
with identical parameters has shown agreement to the fifth 


and sixth decimal places. 


E. SKEW RAYS IN GRIN 

GRIN skew rays are handled analogously to the homogeneous 
case with the same coordinate transformation from grid plane 
to global coordinates being required. It is only after the 
initial directions cosines of the GRIN skew ray K', L', M' 
are found that the differences between GRIN and HIN appear. 
The only exception to this being the use of Equation (113) 
in Snell's law for refraction at the interface. Since the 
GRIN skew rays display curvature in a plane through Oo: Kees 
L', and M‘' are constantly changing until intercept with the 
inside conical surface. Therefore, not only must the plane 
of the skew ray be analytically described, but the final 


Paeseceor K', L', M' must be found. Due to the nature of 
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@arn tays, the procedure for finding the ray intercept with 
ne conical inside surface is different from both the HIN 
case and the iteration procedure employed in the meridian 
plane due to the multiplicity of the geometry encountered. 
The Newton-Raphson iteration routine has been found to be 
ideal for this purpose. 

To begin, the magnitude and direction of the initial 


Zz. must be ascertained. 


radius vector r. from 0. for x 9 


0 
The magnitude is given by 


‘a a 


eZ 


ee cee 6B) eye # zo (145) 


0 0 o} 


Mmerefore the unit vector in the direction of the intercept 


is 
k (146) 


Mid plane of the ray may be fully described by the vector 


normal to the plane. Two vectors, r, and R lie in the plane 


0 


of the ray. R is the unit vector in the initial direction 
of the ray after refraction and described using Gireceion 


cosines as 
R Se US cee ap a ul (147) 


Thus, the plane of the ray may be G@esentbedeoye ene —eross 


product 
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eZ > 
I 
ty > 
x 
wv 


Po Oy (148) 
er 
=e i om jy +N. Ok (149) 
Po Pox Poy Poz 
where 
= Lae ' 
>. Foy Loot (150) 
N = ieee) = Sa (1 Sg 
P oy Oz Ox 
and 
N = eli oak * tS 2 
Pos Ox Oy ( ) 
meeematuons (150), (151), and (152) LOy! roy! and Yoo 


jeer tO the x, vy, and z components of r, in Equation (146). 


0 


Furthermore, the angle Yo between r, and R may be found from 


0 
the dot product as 


x oxk ~ Foyt ~ rooM ] Ci s:3 } 


Now the scalar invariant e may be found. Substituting known 


values into Equation (109): 
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<0) tn 
e = e¢e[art as ] ry sin Vo (154) 
where Ro is unchanged from that found during previous calcu- 
lations for the lens shape by Equation (127). 

With the foregoing groundwork established, the intercept 
of the ray with the inside surface may be calculated. The 
Newton-Raphson iteration scheme requires the calculation of 
tme radius vector from o- to the cone by geometrical methods 
and the radius to the ray by GRIN theory. The difference 
between the two radii is then divided by the difference 
between the derivatives of the two functions. The resulting 
quantity is subtracted from the trial angle, Jae in the plane 
of the ray, to give a new trial angle oy: The process is 
continued until the difference between radii is less than 
Rcclgee The first trial angle is measured to a reference 
HIN intercept as if the material were homogeneous since the 
actual GRIN ray curves only slightly. The coordinates of 
fae lee ie eo 


are derived using the HIN equations as before. To obtain 


this HIN intercept point are designated x 


the first trial angle, the scalar product between a and 


re. 1s used. Here 
1H 


5 5 eZ 


2 
G. = {(x.,, + OB) + ae Zs] ES) 


Dekh 1H 
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ba 





rn >: nia mane ©) > aN Weare A ri 
x 1H ; Seri .. sam 
ae let jt 5 k (156) 
1H art sie 

Tiass., 

6s = cos t[r, r. oe tonto ge ekg) Au eee (sy) 

p Ore ai x Ov thy 2a HZ 

where L ax’ t iHy’ and Tay, are Ene x, VY, and 2 components 
oba Lays respectively. 

my “Eauation (117) Is Psstestieucea for (9-9 9) or yield 
r aS required by the iteration procedure. The geometrical 


radius, sy isete~wwse easily acquired. First, note that the 
equation of the plane of the ray inside the lens is given by 


N x ¥7X) -o ) = 0 CES 3) 


pod 


Wee) * Noo ‘2770 


EOy 
Meeonc!y, the equation of the conical surface is given 

Sy euation (75). The combination of the plane of the ray 
gaa tne con® yield the loci of possible intercept points on 
the inside surface. In Cartesian coordinates, the sum of 


Meuattons (75) and (158) is 


+y] + z([N +z] 


2 , = 
KON  g-% tana] + alr 50z N50x*0 


OY 
=n - = ino, 
N oy’ 0 N 50270 0 ( ) 


@aueeston (159) must be transformed into coordinates r and a 
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in the plane of the ray. The transformation is made possible 


WeamG Bouations (120) through (124). Upon substitution, 
and after solving for Eo Equation (159) becomes 
2 
- | eee: 
r oe ao as x (160) 
g 2 4a 2 
Z 
where 
ms 2 Zz Z 2 
A, = By + Cy AL tan a (161) 
Z 
= L L W, 
B. Boy + CT oz ~ ex + 2A, OB tan oO ClEGZ} 
c =  OBN + ae ano + N co meeter ly 4 
2 pox elope poy~ 0 
+ N Z (163) 
poz 0 
and 
m - “36 = op) + nk" (164) 
id ro 0 
SY 
3 _ ane (165) 
i ro 
525 
C Z + nM! (166) 
ak “ 
0 
Meee, 6 and 4 are found from 
ae ; i no? 
ae = Sin Daf ae bo ( ) 


eal 





anid 
c= eeCcos 2s - COS Vy (168) 


which follow from Equations (123) and (124). The plus sign 
in Equation (160) yields the correct values. Now that r of 
the ray and ry of the surface are known, the derivatives of 
r and ty with respect to 9 at a must be found. It can be 


shown that for the ray 


dir ame =1L 
AG, = Seam. cos[-2.¢90 + sin (A,) | Cres) 
9 2e P 
2? 
where 
Zi Jaa - a 
A, ee ee Eo) 


ac - 4be“/RS 


The derivative of Se is somewhat more complicated. With 


persistence, however, it can be shown that 


2 
dr dB aA B Cred 2 
SS ae _ 2 ee | Se 
a6 me ae 2 Ge) + a! at x [c.] 
2A. 4a, 2 
(171) 
where Ao, Bos and Cc. are given by Equations (161) through 
(168) and 
iB GA C. aA 
z ipe2e M2 2 a ar” ob 
C., = me ad am - Bo Ay iam ae a5! ( ) 
2 2 


oz 








Preresore, the r and Fs to the intercept are found by 


iteration of revised trial values 


(r-r ) 
: = ea 
a5 as a. aE (07) 
ar qe” 


Bi practice the quotient of differences in Equation (173) 

is reduced by a factor of 1.3 to slow convergence and provide 
stability. The number of iterations, however, rarely exceeds 
five. 

Values for r and aS are now transformed into Cartesian 
coordinates Xo, Vir 2, DySsebstiEution into Equations (120) 
through (124). Next, the values of the direction cosines 
K', L', and M' at the intercept are needed. The angle yp 
between the radius vector and the tangent to the ray at 
intercept may first be deduced from the scalar invariant, 


e. Hence 





) (174) 


Where e is known, n. is evaluated at r is the 


2 De Ea 


radius to the intercept as found by the iteration above. 
Three constraints on the direction cosines may be written. 


These are: 


“A 


Tor and the 


ieetne Scalar product of the radius vector, 
Migeekitancous ray direction vector, R, at intercept. 
2) The scalar product of the normal to the plane of the 


ray N50 and R. 


a3 





3) The sum of the squares of the direction cosines must 


Sumto unity. 
Mathematically, the above constraints are written as 


(Xx, + OB) Wee Z 


K' + == L' + == Mm! = cos v (175) 
PI PI PI 


Noo xk + ae + Ns (176) 


I 
© 


> ae +M = 1 a7) 


io oene &, &’', and M', Equations (175), (176) and (177) are 





solved simultaneously. It can be shown that the solution 
leads to 
2 
5 B Ss 
Mis ye t V3-S (178) 
6 4A 6 
6 
C= . - Bm? (12) 
5 5 
Wee 
ei = Li - ae) (180) 
where 
Npoz,2,,2 _ *Npoy “poz 
aoa (le 6) 2 Es 
50x : N 
P pOx 
N502, 2 
+ (P)% + 1 (181) 
pUx 
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2N N N 
- _ pOy  p0z _ Pp 
Be Noo = A, - Jil + (; 2) “1ASBs (182) 
DpUX mux 
nl le me 
Cae il + (E*) JA, - 1 (183) 
pO 


and 


N cos wv 
Ox 
N 50x? 4 Nooy4 


B= “pox 4 ~ "p02" 4 (185) 
pOx 4 p0y 4 

Purthermore, in Equations (184) and (185), Ayr Ba and Cy 

are the coefficients of K', L' and M' in Equation (175). 

Now that the direction cosines of the ray are known at 
the point of intersection with the inside surface, the index 
Gee retraction, No» Ws"comeuted by substituting Cor into 
Equation (113). The remainder of the skew ray trace to the 


imace plane is identical to the homogeneous procedure. 
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IV. LENS PERFORMANCE PARAMETERS 


The function of the seeker lens is to focus electromag- 
netic energy either reflected from or emitted by the target 
onto a detector. Angular displacement of the target with 
respect to the missile body axes as well as target angular 
rate information are both desired outputs from the seeker. 
Mmeerece, Me is net only important just to be able to detect 
the target by focusing energy into a spot on the detector, 
it is equally important that this spot be as small as possi- 
ble to enable the precise position of the spot on the detec- 
tor to be discerned. 

Tne ability of a lens to focus an object to a small spot 
Goes not guarantee the quality of the image. For a FLIR 
Optical system extensive effort is expended to obtain an 
image with minimum aberration [10]. Seeker optics, however, 
are generally non-imaging devices where the pressure of the 
Gifferent aberrations does not detract from the function of 
the seeker as long as a tight image is maintained [ll]. 

Accordingly, the most important parameter by which seeker 
H-ns Perceormance is judged is that of spot size at different 
obliquities. Since the image found at the focal point is not 
necessarily circular nor equally dense, the standard devia- 
Hom in the y and z directions with respect to image centroid 
me Weed to define spot size. 

Mieisetore, once the coordinate pairs of all the skew rays 


have been calculated in the image plane, the first step in 
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the analysis of lens performance is to calculate the image 
centroid [12]. If the spot diagram is composed of N rays, 


the centroid location, a 1s 
N 
1 


Since there is symmetry about the x-y plane, the z coor- 
dinate of the centroid will always be zero. 

The next step in finding the spot size is to find the 
standard deviations of the spot diagram in the y and z 
directions. This is accomplished by summing the squares of 
the differences of the intercept coordinates with respect 
to the centroid and then dividing by the number of rays. 


The standard deviations then, are given by 


N 

Ze, te 2 ps 

tal 3 87) 
i=l 

and 

N 

2 LL = 2 

oe ~ hae ee (188) 


Zhe Spot size, Sh is now defined by 


m= “los = o- (189) 


my 





Of further interest in appraising lens performance is 
mee €Rengy Censity of the image as a function of radius from 
eae centroid. In nondimensional fom, this is simply the 
number of rays in the spot diagram within a succession of 
Circles of increasing size overlaid about the centroid. Here, 
each ray 1S assumed to carry a unit amount of radiant energy. 
This type of plot facilitates the comparison of different 
lens designs by detailing the distribution of energy within 
each image. Clearly, it is desirable to have as much energy 
as possible concentrated very close to the image centroid. 
Between two lenses with equal spot sizes, the preferred lens 
has more energy concentrated within a smaller radius. 

Bvery Gay, hewever, does not deliver an equal amount of 
Puemeay to the focal plane. It is prudent, therefore, to 
include as a performance parameter the average ray intensity. 
iaaim, cor N rays, the intensity I of each ray 1s summed and 


Hesmalized by N to yield 


a 
iE Si ) ae (190) 


Finally, each lens design is checked for "failed rays". 
The reader will recall that these rays fail to intercept the 
image plane due to total internal reflection, total external 
reflection, or failure to intercept the inside surface within 
the bounds of the lens. Hence, a lens design with fewer 


Meauiledszays" or no "failed rays” at all is a preferred lens. 
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Although there are many other performance criteria by 
which lenses are compared, the foregoing parameters are more 
than sufficient to judge the merit of preliminary seeker 
lens designs. It should be noted, however, that notwithstand- 
ing the fact that the image centroid and standard deviations 
were used as stepping stones to obtain image spot size, they 
have significant meaning of their own. The standard devia- 
tions a and e. inform the lens designer as to the horizontal 
and vertical spread of the image. Image centroid location, 
a! at increasing obligquities is of obvious importance since 
excessive displacement will cause the image to miss the 
detector entirely and would dictate the necessity for a 
second lens element to dampen the movement. Furthermore, 
Since Line of Sight (LOS) measurement accuracy to the target 
is highly dependent upon the linearity of Y, as a LNneGe Lon, 7.O. 
the lens tilt angle, lens designs which exhibit a greater 
degree of such linear behavior are the preferred designs. 

To sWmmarize, it is sufficient to note that although 
SeOct Size is the most important of the performance parameters, 
every other parameter has a significant impact on the 


mertormance of a particular lens. 
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V. RESULTS FOR THE HOMOGENEOUS LENS 


The performance of the homogeneous lens is presented 
primarily as a comparison with which to compare the performance 
@2 ehe GRIN lens. Here, the relationship of spot size to 
increasing lens obliquity and lens thickness are presented 
as well as the image centroid movement as a function of 
obliquity. Additionally, the reader is introduced to the 
four basic computer plots used to display the results: lens 
shape, object plane with superimposed skew ray grid, image 
plane or spot diagram, and energy density. These plots were 
generated on the VERSATEC Plotter using arrays of data points 
produced by program GISL on the IBM 3033 mainframe computer. 

To begin, Figure 1l shows the homogeneous lens shape. 

The first of the four basic plots, the lens shape plot, pre- 
sents the lens side view in the meridian plane. The outer 
Surface (curved) and the inner surface (cone) are constructed 
by connecting the points x, (d) , y, (A) and X(T) py¥5(I), 
respectively, by straight lines. At the apex of the lens 
about the lens axis iS a trapezoidal region which represents 
the opacue nose area. In the legend are listed the lens 
design parameters and the significant calculated dimensions 
of the lens. Since all linear dimensions are implicitly 
normalized with respect to the maximum inside radius, R, 

tne lens may be scaled up or down by multiplying each dimen- 
Swen by a factor of RA /R. Beginning at the top, the 


Cw 


parameters listed in the legend are explained in Table l. 
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Figure 11. Homogeneous Lens Shape fOr No = iD 
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PARAMETER 


F/ # 
E 


OS YMB 


ALPHA 


B 


Serr 


% Index 


% Inside 


% Outside 


TABLE l 


Explanation of Lens Shape Plot Legend 


TYPE OF PARAMETER 


DESIGN 


DESIGN 


DES IGN 


DES IGN 


DESIGN 


DES IGN 


ESE Gil 


DESIGN 


DESIGN 


DESIGN 


DESIGN 


DESIGN 


CALCULATED 


DESIGN 


CALCULATED 


CALCULATED 


OZ 


MEANING 


Fo number. F/# = F/2R 


* 
Focal length from B. 


OB. Line segment O_ to B 


(positive to left) 7 
Immaterial in HIN 


* 
Maximum radius of cone 


* 
a--cone half-angle, radians 


* 
Edge thickness 


Incident ray offset angle 
(design) *, radians 


Number of iterations. 
I¢+l = number of rays 


Free stream index of 
refraction 


Pndex<.Or set racktoneo dt 
interior lens cavity 


a in n(x) = Va+b(x°/R, 


gradient refractive - 
index function. N, = Va 
me HN. 


joeeene| no(r). Zero in HIN 


b effective = b/RE 


Percent change in n. (zx) 


from r = 0 to r = Ro 


(73) 


fom 


Percent change inn 
along inside surfacé 
lens axis to edge 


Percent change in n.(r) 
along outside surface 
from opaque region to edge 





rable 1 (Continued) 


PARAMETER ZYPE OF PARAMETER 


% Across CALCULATED 
STATION A CALCULATED 
GAMMA CALCULATED 
TSURFL CALCULATED 


k 
MEANING 


Percent change in n.(r) 
across lens from lens 
axis to outside surface 
at the thickest point 


xX-coordinate at nose of 
lens 


y--nose half-angle of 
Opaque nose region, 
radians 


Total outside surface 
length from Station A 
to the edge 


“Refer to premee 2 for clarification 
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The HIN lens shape has a convex outer surface with maxi- 
mum thickness on axis of almost ten times the edge thickness. 
Although the lens has a good aerodynamic shape resembling an 
Ogive, the outer surface is not a circular arc nor can a 
Single analytical function be fitted to the array of points 
Gescribing the surface. Note that the nose half-angle, y, 
is almost identical to the cone half-angle, a. 

All lenses have been designed with a cone half-angle of 
45° which is approximately the maximum angle for which aero- 
dynamically efficient lens shapes may be designed, consider- 
ing a free stream Mach Number not to exceed three. Without 
exception, overall lens performance is more severely degraded 
as angle a is reduced. 

Table 2 explains the legend of Figure 12, which is the 
second basic plot. Here the lens is depicted as seen from 
the skew ray grid plane. The lens tilt angle, a5? causes 
the equally spaced (in J) circles descriging the surface of the 
lens to appear as ellipses. In Figure 12, the grid spacing has 
ween reduced from 0.1, which is normally used, to 0.3 to 
allow identification of individual rays for correlation with 
the image plane spot diagram. Although the small number of 
rays used is not sufficient to given an accurate definition 
of spot size, the number is sufficient to describe where 
rays in the object plane are being focused in the image plane 
by the HIN lens. The skew rays in Figure 12 have been numbered 
in the order in which they were processed. Actually, only 


fms lthrough 19%were actually traced; 20 through 33 are 
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TABLE 2 


Explanation of OBJECT PLANE Plot Legend 


PARAMETER TYPE OF PARAMETER MEANING 

ALFAP Analysis Or) LeNerecd le vag le, 
rSdians 

Rays Analysis Total number of rays 
processed 

Failed Rays Analysis Total number of rays 
igs eb iLabgye: to pass Ennough 
the lens 


* 
betbed Rays are indicated on the plot by a diamond 
Superimposed on the grid location of the ray. 


mirror image skew rays. Ray 30 corresponds to ray 14, for 
example, Here, the lens has been tilted by 0.4 radians, or 
22.9 degrees, and 33 rays have been processed of which none 
have failed to intercept the image plane. The staircase 
pattern has been added in this case in order to show the 
resulting distortion present in the image plane (Figure 13). 
The Spot Diagram in Figure 13 1s an example of the third 
basic computer plot; see Table 3. Unlike most Spot Diagrams, 
this example has the individual rays numbered for comparison with 
Figure 12; also the resulting distorted staircase pattern 
is sketched. By cross-referencing individual rays between 
Brmeeres 12 and 13, it is possible to recognize where certain 
areas of the lens are focusing rays in the image plane. Rays 
1, 4, 20, and 22 about the opaque nose region form a coma 


tail which contributes most of the image spread. Rays 33, 
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= ie 
33 19 
1 x 18 * 
. 
6 AO Mm 14 24 
we § pers? 27410 4 \a3, #21 
7 
25% mm *" C8 
26 3 9 
+e 
ae «3 — = 
; FAP & 0.400000 
SPOKASIAE =0.4987394 
a2 Y CENTAWQ = -1.5702500 


SIGMA _ =". 2029376 
S GMA 2 =O. 0456086 
RAY COUNT =33 





Figure 13. Spot Diagram for HIN Lens Design Shown 
in Figure ll 
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TABLE 3 


Explanation of Spot Diagram Legend 


PARAMETER TYPE OF PARAMETER MEANING 

ALFAP Analysis tn plencwetltivenale, 
radians 

SPOT SIZE Performance O., see Equation (189) 

Y-CENTROID Performance Yor y-coordinate of 


image centroid; see 
Equation (186) 


SIGMA Y Performance Ov, y-Standard deviation; 
Seesrquation (133) 


SIGMA Z Pestommance Og, z-standard deviation; 
see Equation (187) 


RAY COUNT Performance Number of rays striking 
the image plane 


18, and 19 from the bottom portion of the homogeneous lens 
are imaged at the top and are widely separated from the core 
@i the image. In general, the regions of the lens which 
have been found to contribute the bulk of the widely spaced 
rays are the immediate nose region and the lower portion of 
the lens. 

Furthermore, the distortion present in the image of the 
staircase pattern clearly shows that regions closest to the 
nose yield the greatest distortion. Horizontal lines are 
Switched end for end and tilted approximately 45 degrees. 

Tee Upper portion of the lens performs the best. Rays 
See, et, and 15 an the meridian plane are focused on the 


You “amor a faarly tight region; rays 10, 14, 17, 27, 30, 
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and 32 about the upper periphery in the object plane are all 
focused within the image core. 

Table 3 explains the legend of the spot diagram. It is 
seen that the standard deviation in the y-direction is 
approximately five times that of the z-direction. This 
elongation of the image is not readily evident in the spot 
diagram since the ordinate and abscissa have not been plotted 
with equal increments. This results from the great disparity 
between Spot Diagrams of the various lens designs studied. 

It 1S important, therefore, for the reader to take careful 
note of the relative sizes of the Yim and Zrm axes. 

Of primary importance is the spot size in Figure 13. The 
value for spot size is adversely affected by the poorest 
performing regions of the lens. Were it not for these 
errant rays, the spot size would be considerably smaller. 

The image intensity pattern is benefitted, however, by the 
fact that the rays spread the farthest from the centroid 
contribute significantly less energy per ray than those 

being focused in the core of the image. Table 4 lists the 
relative intensities of the primary skew rays plotted in 
Figure 13. It is seen that the high intercept angles experi- 
enced by the rays closest to the bottom of the lens and, to 

a lesser extent, those near the nose, result in higher 
reflectivity and lower transmission through the lens. 

In the legend of Figure 14, the average of ray intensity 
is given. The Spot Diagram Energy Density distribution may 


be seen at a glance. The fraction of energy (number of rays) 


ee 


TABLE 4 


owew Ray Intensities of HIN Lens 


a) Nwmberical Order 


Ray Dmeensacy Ray ineenst by 
1 0.746 11 Os StS) 
2 Gee yatens) ip OS thes 
x) Ono 21 1-3 OpeoD 
4 Oi, Gore 14 Omccn 
5 OR847 Ibs OF oie 
6 O7-o I NG OS) ILS) 
au Orga 0 1 0.910 
8 Om 7 18 0.640 
9 0 2904 as 0.576 

Ke 02870 


b) In order of Descending Intensity 


Intensity Ray (Ss) ingens. ty Ray (s) 
O.921 0.887 14 
0.920 r, Oe eter 5 
e719 ies 0.870 ie, 
@.915 Id es 0.811 
OR 7 8 O27 59 
02916 16 0.746 dt 
O29 10 ie oe ly 0.640 18 
0.904 £ O57 6 ths, 

O #69) 2 4 
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Shown as a function of radius from the Spot Diagram centroid 
is normalized with respect to the total number of rays 
traces (Rays, in Figure 12) regardless of whether all of 

the rays successfully intercepted the image plane. Thus 

the fourth basic plot may be used in conjunction with the 
Spot Diagram to further define the image concentration with 
respect to the centroid. 

The response of the HIN lens to increasing tilt angles 
1s given in Figures p-1 through D-23 where the index of refrac- 
tion of the lens has been set at 1.5. Whereas the spot size 
at 8S = 0.0 is very small, that at a = 0.7 radians is quite 
large at 67% of the lens radius. Figures 15 and 16 summarize 
spot size and centroid locations for the lens. Spot size 
growth, although somewhat irregular at the higher angles, 
1s pronounced. Furthermore, the centroid movement is seen 
to be an approximately linear function of o ee 2 1 
a = 0.3 and easily exceeds the radial dimension of the lens 
mee igner tilt angles. 

Figure 17 shows that the HIN lens may be slightly improved 
by increasing the edge thickness, T, by a small amount. The 
lens may be otherwise tuned to improve performance at certain 
tilt angles by designing the lens with U slightly greater 
than zero. These performance improvements are practically 
insignificant, however, and neither lens tuning by the param 
eter U nor T produce improvement across the spectrum of tilt 


angles. Instead, an improvement at one ore usually has resulted 


in a degradation at others. 
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meee 17. HIN’ Lens Spot Size versus Edge Thickness 
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Finally, the homogeneous lens is shown with index of 
reteaction of three at the intermediate tilt angle of 0.3 
radians in Figures D-24 through D-27. The lens shape required 
to accommodate the higher index of refraction is seen to be 
thinner and displays less outside surface curvature than 
the HIN lens with N, elope ODOrmSUZ7C US .SIighisiecantly re—- 
duced. Since infrared lenses such as germanium generally 


have fairly high refractive indices, improved lens per- 


formance at these higher values is encouraging. 
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VI. GRIN LENS RESULTS 


The performance of the GRIN lens is similar in many 
respects to the homogeneous lens. The relationship between 
characteristic regions of the GRIN lens, such as the lower 
lens portion, and where these regions image bundles of rays 
is identical to the HIN lens as depicted in Figures il and 
12. The growth of spot size with increasing or and the 
respective movement of image centroid typical of the homogene- 
ous lens is clearly displayed by the GRIN lens as well. 

The measure of the superiority of the GRIN lens, there- 
fore, lies in the successful correction or improvement of 
the deficiencies seen in the HIN lens. Here, the reduction 
of spot size is of primary concern. 

As a modest example of the ability of the gradient refrac- 
tive index to reduce spot size, the GRIN lens design shown 
in Figure 18 is examined. This lens is very Similar in shape 
to the homogeneous lens with No = 1.5. Note, however, that 
unlike the HIN lens, the nose half-angle, y, 1s slightly 
larger than that of the cone, indicating more outside surface 
curvature. The object plane for this lens, at os = 0.4, is 
Shown in Figure 19. Here, the error in focusing rays present 
in the image plane has been superimposed over the grid plane. 
You coordinate error contours (with respect to the centroid) 


are shown in the left half of the plane and Zou contours in 


the right. These contours vividly show that the largest 
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Figure 18. Example GRIN Lens Design with 10% 
Negative Gradient at OB = 1.0 
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errors in both y and z stem from the nose and bottom of the 
lens. The best performance is contributed by the upper 
central region. The reader will note that negative errors 
imeZ ocemr below the Yo = 0 plane of the lens and that nega- 
tive errors in y straddle the Soe O plane. 

The corresponding intensity contours are displayed in 
Figure 20. Of significance here is that regions which perform 
relatively well in y and z error also perform well in trans- 
mitting energy. Regions of the lens which are characterized 
by relatively high angles of incidence, therefore, perform 
the poorest and have the most to gain from better combinations 
of gradient index and Og 

Mie "spot Diagram, Figure 21, displaysjan image pattern 
wiremally identical with the HIN lens in Figure ll. The 
GRIN spot size, Os however, is slightly smaller by approxi- 
mately 2% and both oy, and c¢, are correspondingly smaller 
which indicates superior performance by the GRIN version. 
Image centroid for the GRIN is displaced further than that 
for HIN. The average intensity transmitted by the GRIN lens 
compares favorably at 0.88 as opposed to the average inten- 
sity of the HIN at 0.85. Further evidence of superiority 
Gpeeie "Gnls Pens 1s seen in the form of a steeper slope to 
the Encircled Energy Plot (Figure 22) than that of the HIN 
tens in Figure 13. 

Although the example GRIN design displays marginal 


Superiority, some GRIN designs do not. [In particular, it 
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Figure 22. Example GRIN Encircled Energy Plot 
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has been found that as the center of symmetry of the index 
function, Oc, is moved further away from the interior of the 
GRIN lens, performance deteriorates. 

Two series of GRIN lens designs have been examined at 
the intermediate tilt angle of 0.3 radians in order to define 
the spectrum of lens performance. Both design series explore 
the GRIN lens for 0. inside, outside, and far outside the 
interior of the lens. The first series, Figures E-l through 
E-92, are GRIN designs in the "low range" of refractive index 
with the index parameter "a" set at 2.25. This first set is 
compared to the HIN lens with = va = 185. “Ehe second, 
or “high range” series is for a = 9.0, as compared to the 
Him at No = 3.0, and are displayed in Figures F-l to F-128. 
Both the high and low range lenses exploit gradient changes 
or 5, 10, 25, and 50 percent both positive and negative, 
where possible. A negative 50% gradient in the low range 
is, of course, not possible since an index of below 1.0 
would result. The center of symmetry was not located at 
x > 0 due to the aforementioned singularity encountered at 
very small angles. 

The lens shapes of the GRIN designs shown vary widely. 
Tne outside surface may be either convex or concave and ex- 
hibits a thinner profile at higher values of refractive 
m@eex. Although all of the resultant GRIN lens shapes are 


Superior aerodynamically to the hemispherical seeker lenses 


currently in use, the convex version of the lens has more 
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obvious applications for a Sidewinder type tactical missile, 
whereas the concave lens shape is more applicable to the 
diffuser of a ramjet with nose inlet. 

Object plane diagrams show that for lens designs exhibit- 
ing smaller spot size, fewer or none of the skew rays failed. 
This fact underscores the success of GRIN in controlling image 
@eeerioration contributed by the nose region of the lenses. 
Although the lower lens and nose regions still create the 
largest of the image spread, the extent is reduced. 

Figures 23 and 24 are contour diagrams summarizing 
spot size performance of GRIN lens designs in the low and 
high range respectively. From these diagrams, it can be seen 
immeciately that the best performance is obtained from lenses 
having the center of symmetry inside the lens. From this 
fact, it may be deduced that large changes in refractive 
index are desired along the surface of the lens rather than 
across. Except for isolated regions at -5%, increasing 
Cradsenmts Daro@uce smialler spot sizes. The positive gradient 
in the high range, however, has an almost constant spot size 
from +5% to +50% where an improvement of only 0.6% is seen. 
In both the low and high ranges the positive gradient at 
50% and at OB = 0.05 has proved to be the best performing 
combination. Again, as in the HIN lens the high range 
exhibits the smaller spot size. 

The best GRIN lens is now further examined over the range 


Sieebliguicies up to 0.8 radians. Furthermore, since it has 
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been determined that placing the focal point at x = 2.0, 
improves spot size performance at lower values of Oo another 
series of plots is given for F/l. With the focal POLnt at 

x = 1.5, severe degradation results; for F = 3.0 Only modest 
differences with F/2 are evident. Appendices G and H show 
the pexptormaence of the "best" GRIN lens at F/1 and F/2. 

Here it may be fully recognized that the resulting be- 
havior of a "good" gradient refractive index seeker lens is 
far superior to that of the HIN lens for Os < 0.6 radians; 
see Table 5. The reader will note that there are no failed 
rays even at 0.8 radians (45.8 degrees). Although centroid 
movement at F/2 is slightly greater than the HIN lens, by 
changing to F/1l, this can be corrected; see Figure 25. The 
F/l version, moreover, displays a range of spot sizes signi- 
ficantly smaller than that for the F/2 lens below O05 = 0.54 
radians or the HIN lens below 0.68 radians; see Figure 26. 
That the lens displays more desirable performance with a 
Shorter focal length is a surprising but highly desired re- 
sult. Since seeker optics systems are volumn limited, any 
reduction in the lens focal length is beneficial to the final 
Gesign and vackaging requirements. 

It is significant for the reader to note that centroid 
movement beyond a value of +1.0 exceeds the physical dimen- 
Sions of the lens radius. If the interior of the lens mount 
is cf similar dimension, it follows that the requirement for 
Yoentr to be less tiemeeai Ul restricts the F/] lens to approxi- 
@aeely 0.65 radians of tilt whereas the F/2 lens is further 


restricted to less than 0.3 radians. The superior performance 
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TABEE 5 


Comparison of HIN and GRIN Lens Designs 


* 
Performance Parameter or Feature HIN GRIN 


Shape of Outer Surface Convex Concave 


Fraction of Failed Rays 
at = = 0.7 radians 0.143 0.0 


Linearity of Centroid 
Motion; a, for Ay/y 0.4 radians 0.8 radians 
deviation of 103% 


Spee Saze for “ss 


- =) @.1 waditans 0.190 0.029 
cee = 0.2 radians 0.247 0.063 
ag =s0.3 Padrtans 0.442 Oe aleie2 
- =— 0 2 radians 0.466 02201 
= =—eQe5 radwans 0.637 On 33 
a = 0.6 radaans GG. 2 0.490 
a = 0.7 radians 0.667 0.688 
= = 0.8 radians N/A 029 47 
Seoe eacius for 30% 
of energy @a_ = 0.4 Does not 0.28 
radians P attain 803% 
waetnin 0.50 
Average relative intensity 
of skew rays 34 Cag G4 0.82 Oa 


radians = 


* ° « 
"Best" GRIN lens with 50% positive gradient, 0. OO) Sie 
ia the F/1 configuration. 
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of the F/2 lens above 0.54 radians 1s, therefore, unusable. 
Attempts to improve upon this performance by increasing the 
edge thickness, T, or by slightly adjusting the parameter U, 
as with the homogeneous lens only produced uniformly degraded 
performance in every respect. 

Although the "best" GRIN lens resulted from a positive 
90% gradient change in refractive index, it has been noted 
Ghiae this cemftiguration was only slightly better than the 
same lens with a positive 5% change. Furthermore, since 
refractive index gradients of five percent or better have 
already been produced, it is entirely feasible that if the 
precise parabolic change could be controlled, this lens could 
be produced today. 

Despite the obvious success of GRIN in controlling spot 
size growth and image centroid movement, a penalty in the 
form of reduced ray intensity has been paid. That increasing 
spot size performance is tempered by a loss of intensity may 
be seen by comparing "good" GRIN encircled energy plots with 
that of the HIN lens. Note that the HIN lens with No = 3.0 
also loses ~ 50% to intensity; see Figure 44. This loss in 
intensity is partially offset, however, by the GRIN lens 
with an increased number of rays transmitted and by the re- 
duced spot size. 

One drawback to the use of the GRIN lens as a self 
sufficient, single element lens is still that of the relative 


sizes of image and detector. Hence, even though GRIN has 


lig 2 





ile 


Significantly reduced spot size below os 0.6, that size is 
still significant at tilt angles above 0.3 radians. In order 
to use the lens without a secondary focusing element or 


Mirror arrangement, a large composite sensor array would be 


required. 


IV eyes: 





Vite eCONCLUS LON 


Gradient refractive index materials may be employed to 
design a pointed seeker lens which exhibits optical per- 
formance far superior to that obtainable with conventional 
homogeneous optical material. A fifty percent, positive, 
parabolic gradient index with center of symmetry interior 
to the lens was found to yield the best performance although 
a five percent version of the same lens was only very slightly 
inferior; this lens may possibly be fabricated today. 

For objects off-axis by more than 0.3 radians (~17.2 
degrees) a secondary lens element may be required unless a 
large scale multiple element sensor array is employed. With 
Sue weanmays»uob jects off-axis by more than 0.65 radians (37.2 
degrees) may require Cassegrainian or other mirror elements 


to compensate image movement. 
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VIII. RECOMMENDATIONS FOR FUTURE WORK 


This thesis investigated a spherically symmetric GRIN 


seeker lens with inside conical surface. Future studies 


should investigate: 


a) 


b) 


Cc) 


dq) 


e) 


£) 


The elimination of the singularity problem by 
modification of the theoretical equations for GRIN 
using numerical techniques. See Reference 13. 

The optimization of the GRIN lens with both inside 
and outside curvature. 

The effect on lens performance due to an attached 
shock wave. 

The effect on lens performance when the object is 
no longer in the far field and wave front curvature 
must be taken into account. 

The performance of the lens using wavelengths corres- 
ponding to atmospheric windows. 

The effect of a radially symmetric gradient index 
on lens performance. 

The feasibility of adding an anamorphic gradient- 
index lens located on the missile body to increase 
@fft-aems tracking/acquisition Capability to 90° and 


beyond. See Figure 1 of Reference 14. 
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APPENDIX A 


COMPUTER LOGIC FLOW DIAGRAM FOR PROGRAM GISL 


This appendix together with Appendices B and C describe 
the FORTRAN program GISL. GISL may be used to design either 
a homogeneous (B = 0), or a GRIN lens to the user's specifi- 
cations by changing the design parameter where indicated by 
comments in the input section of the program. Additional 
software required to run GISL and plot the results are not 
described since these system-dependent procedures do not 
apply elsewhere. 

The flow diagrams which follow provide the reader with 
the information necessary to follow the fundamental logic 
of the main program and subroutines of GISL. Not shown is 
Subroutine DIRECT which has been derived from Subroutine 
DEeeeeas described by Amichai [5] with only slight modifi- 


cations. 
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INDEX, 
RADIUS (3) 


x1 (1) ,¥1(1) se 
DYDXN (1) Ces 
SECTION 





DY DXT (1) 


eS 
COMP 
CESS Giles 





COMPUTE 
DYDXN(J) , 
By bx (a) 








X1(J) = X1P 
Winks} = YIP 


X1(J) = X1H 
Y1(J) = Y1H 












ou 





INITIALIZE INCREMENT 
RAYY (G) RAY Z (G) 
RAYZ (G) 





SEARCH FOR 
Y1 (J) >RAY (G) 





PRINT CALCULATE 
500 R 





RAYY (G) ; greiner 
oe!” ~~) payz(G) =0 





* 
Here, the prime indicates Y-coordinate as transformed into the grid plane. 


3 








GRA t yan 


Ea 
INDEX, 


DIRECT 


Y1M,Z1M 





MIRROR 
SKEW RAYS 
MATCHED 

ELLIPSES 
YCENTR 
CALCULATE 


VERSATEC OR 
PRINTER 


NO 


PLOTS 


END 
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CALCULATE 
RPER, RPAR 


CALCULATE 
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CALCULATE 
ene Celis 
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Disk 


2 


AP 
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CALCULATE 
A,B,C 


CALCULATE 
A2gB2, C2 
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APPENDIX B 


LISTING OF FORTRAN PROGRAM GISL 
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gemed fume FE PE peed Pag frees PP PED femang Pred feed gemng Deed penny aE Geet pean a Pred permng PHME fey PE med freed Fd remy FE Peng prencey Od fad Fd Deed pemng frm Fd femny PE fem FE Pend geeng PT Ped peewee Pt 
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SOO STOCOCOOCSCOO COCO SCOSCOCOCOCOCOCOODOOCOCOOOSCOOCOOCOCOOOo 
NOPD DOO SANM|M PMO OM OHO SOMO DOMH DHOAAIMGINO™ OF OMOQUMTPNOM OM 
DLAWUTVLIAR LAWN SO ODDO OO ODODOE EE EEE EEE OTCMODCDOSODDOAKRANANOTNDAM 
Se ae ee I OS OOO OO OOOO DOO COOSC ODO COO OOCOCOOCOOSDO Oooo o 
re a es a sO OOOO OO OC OCC COO OOCOCS COO O OOO OOO OOOO CCOOoO 
NMNNNYViNWNIDN HYMN MNNNYVIN MNNNMN NUM NY NNWUINNNMNNNNNMYLNYMMMNNMWN 
81565508 0555.5 5 COON COLON CD PRED EH CHED ENED COCHED COCO EDO 


SQOOQOOO90S eee = OO OW OOOO 0Oree EHO OOO UICOOOOO 


OTHERWISE 


USHAPE” TO 1. 


SET 


ENS SHAPE CATA PRINTED 


WO OUO YVOO OVOUO O LO 


TA PRINTED SET 


RAY DA 
(INTEGER ) 


IMAGE SKEW 
TO 0. 


ELSE SET "SeuRRY" 


TO 1. 


TO HAVE SKEW RAY AND MIRROR 
wR AY. 
SQURAY=1 


OTtERWISE 


le 


T9 


Shilpa? DATA (PRINTED SET “YELLIPS" 


SE 
Oe 


INPUT KNOWNS FOR GLY LENS SHAPE 


LENGTH: 


ANGLE 


2 wwe 


ae NE HALE 
FFSET ANGLE: 


FACE 
5396 
AY 0 


R 
8 
R 
0 


0 
N 
O 


(MUST BE AN EVEN INTEGER GT 750) 


ITERATIONS 
UM RADIUS OF INSIDE 
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Figure D-4. Spot Diagram Corresponding to Figure D-3 
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Figure D-10. Spot Diagram of HIN Lens at os = 0.3 
Radians, N. = 125 
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Figure D-14. Encircled Energy Plot at 0 = 0.4 
Radians, N, = 1.5 
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Figure D-16. Spot Diagram of HIN Lens at os = 0.5 
Radians, N, = 1.5 
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